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I. INTRODUCTION 

Currently, primordial inflation is supposed to be the leading candidate to provide 
the necessary conditions for the successful big bang cosmology The simplest model of 
inflation driven by only a single inflaton field is consistent with most recent observations jsj. 
It is however expected that, in the context of theories beyond the standard model of particle 
physics e.g. supersymmetry, there is a number of multiple scalar fields which may contribute 



to the inflationary dynamics j4|. Furthermore, we may be able to observationally detect 
deviations from the predictions of single field models in the near future and to discuss 
interesting phenomenology, such as isocurvature perturbations and non-Gaussianity. 



Hybrid inflation ^| is an interesting realization with two field contents, the usual inflaton 
field (j) which drives slow-roll inflation and the waterfall field x which terminates inflation by 
triggering an instability, a "waterfall" phase transition. Previously, it has been assumed that 
X becomes momentarily massless only at the time of waterfall and very heavy otherwise, 
and thus does not contribute to the curvature perturbation TZc on large scales: only the 
quantum fluctuations of contributes to TZc and we can follow the well-known calculations 
of single field case, with the energy density of the universe being dominated by a non-zero 
vacuum energy. 



This naive picture has been receiving a renewed interest with the common qualitative 
results that the power spectrum of the curvature perturbation induced by the waterfall field 
is very blue and extremely small on large scales^. However, quantitatively it is not clear if 







the 6N formalism, which takes account of 



they all agree or not. In particular, in Ref. 
fluctuations only on super-horizon scales by construction, was employed to derive the power 
spectrum, but the approach there was not quantitative enough and hence the dependence 
on the model parameters was not explicitly presented. 

In this note, we provide another complementary view. We adopt a few reasonable as- 
sumptions on the model parameters and solve the mode functions of x in terms of the number 
of e-folds analytically. Then using the SN formalism 10l| we calculate the corresponding TZc 
induced by x explicitly. 

The result is consistent with the above references, i.e. the contribution of x to the 



^ For early atterapts, see e.g. Ref. 
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large scale curvature perturbation is totally negligible. We also clarify the model parameter 
dependence on the spectrum of the curvature perturbation. Furthermore, we calculate the 
corresponding bispectrum, which shows its peak at the equilateral limit. We also compute 
explicitly the one-point probability distribution function which clearly shows the highly 
non-Gaussian nature of the curvature perturbation. 

The outline of this note is as follows. In Section |Tll we find the mode function solution of 
the waterfall field x valid both on super-horizon and sub-horizon scales. In Section llllt we 
calculate the corresponding curvature perturbation IZc induced by x using the 5N formalism. 
In Section HVl we present the power spectrum and bispectrum of IZc- In Section IVj we show 
the explicit form of the one-point probability distribution function of TZc and discuss relates 
issues. We conclude in Section |VT1 In Appendices, we discuss some technical details. In 
Appendix \^ to check the consistency of the 5N formalism with the standard perturbation 
theory, we give an estimation of the curvature perturbation by using the linear perturbation 
equation for TZc- We find a good agreement with our result based on the 5N formalism. In 
Appendix [B] we reconsider the splitting of the super- and sub-horizon modes and compute 
the average over the horizon scales. The results agree with the formulae we use in the main 
text. 

II. MODE FUNCTION SOLUTION OF WATERFALL FIELD 

Before we begin explicit computations, first of all we make the physical picture clear. Our 
purpose is to calculate the contribution of the waterfall field x to the curvature perturbation 
TZc- This is only possible when x becomes dynamically relevant. While x is well anchored 
at its minimum during the phase of slow-roll inflation and hence does not participate in the 
inflationary dynamics, it controls the physical processes from the moment of waterfall till 
the end of inflation. Thus, in the context of the 5N formalism, if we can find the evolution of 
X during this phase as a function of the number of e-folds A^, it amounts to finding TZc by the 
geometrical identity TZc = ^N- Therefore, our aim in this section is to calculate x = x{^) 
starting from the moment of waterfall. We will directly use this result to calculate TZc in 
the next section. 
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We consider the potential of the two fields, the infiaton and the waterfall field as 

^(0, ^) = \{^-^")' + l"^"^" + ■ (2-1) 

We note that during the most period of infiation of our interest, it is assumed that the 
vacuum energy Vq = M^/(4A) dominates so that the Hubble parameter is effectively a 
constant, H = Hq. This is a good approximation even after the waterfall phase transition 
until the last moment of inflation. The slow-roll and the waterfall conditions are 

2 

m 

772 «1' (2-2) 

^0 



M 



Ho 



(2.3) 



respectively. 

The equations of motion are given by 



+ 3H^ + [w? + g^y^) =0 , (2.4) 
X + 3/7x - ^ V^x + (-M^ + g^<^^ + Ax') X =0 , (2.5) 

where the spatial gradient term for is neglected as usual. Note that before waterfall, 
0^ > 0^ = M^/(7^, X is well anchored at its minimum x = so it is itself the same as its 
fluctuation, x = ^X- Thus we may regard (12.51) as the equation for (5x', which arises from 
the vacuum fluctuations. Then after the waterfall transition, 5x becomes unstable and 5x' 
starts to grow rapidly, and inflation ends when the inflaton starts to roll fast, which happens 
when the term g'^5x^ exceeds m? in (12. 4p . Here we adopt the mean fleld approximation, i.e. 
we replace g'^Sx^ by its expectation value g^i^x^), which should be valid for the motion of 
the homogeneous inflaton fleld 0. We also assume that the nonlinear term XSx^ in (12.51) can 
be neglected until the end of inflation. That is, we assume 

M'»4"^'>A(5x'). (2.6) 
9 

At the end of calculation, we must check if this condition is satisfled for the range of the 
parameters of our interest. 



Note that during inflation Spc/, ^ 5(j> while Sp^ ^ Sx^, and thus the metric fluctuations are relatively 
second order with respect to Sx and does not appear in the equation of motion for 6x- This situation is 
closely analogous to the case of false vacuum inflation The correlation functions also show similar 
momentum dependence to those produced during false vacuum inflation 

HQ 
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We can rewrite (12 ■4p and (12. 5p in a more convenient form by using the number of e-folds 
as the time variable, dN = Hdt. Denoting the derivative with respect to by a prime, we 
write 
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0" + 30' +\jj2+ 3 -^) ^ =0 ' (2.7) 

^X" + 35x' - ^V^^x + + + 5x =0 . (2.8) 

Let Nc be the time at which the waterfall transition occurs, (f){Nc) = 0c = M/g. Before 
waterfall, since 6x is very massive, g'^(f)'^ ^ H^, it is dominated by the standard vacuum 
fluctuations and the bare expectation value (Sx^) is ultraviolet divergent. Here we regularize 
it so that it vanishes before waterfall, (Sx^) = at N < Nc. Then (12.71) is easily solved to 
give 

= 0ce-^", (2.9) 

where n = — A^^c is the number of e-folds measured relative to the time of the waterfall 
transition, and we have introduced the parameter r by 



3 /9 m? 

We note that we can write the scale factor a and the conformal time rj = —l/{aH) using n 
as 

a=ace" = ^e", (2.11) 
-no 

-n -n 

^ = -V = -l-' (2.12) 
respectively, where Oc = a{Nc) and = a^HQ. 

Inserting the background solution (12. 9 P for into (12. 8p and neglecting the nonlinear term 
in accordance with the assumption (12. 6p . we obtain the equation for 5x in the Fourier space. 



+ 35x'fc + 



kt 



5Xfc = 0. (2.13) 



A. High frequency limit k/a — > oo: WKB solution 

In the high frequency limit, we can solve (12.131) in terms of the WKB approximation. In 
this limit the proper asymptotic behavior of the positive frequency function is given by 

Sxk = ^ °, exp -2f / dne-" . (2.14) 



2ka .Jl^^kjk^ \ k 
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The WKB solution that has this asymptotic behavior is readily obtained as 



/2M 



TTiexp 



dne 



^ 1 + /3e2- 



where for convenience we have defined (3 by 



/3 = /3(e-"™ - 1) 



The above WKB solution is valid for any k at sufficiently early times, —n ^ 1. 



B. Low frequency limit k/a ^ 0: Hankel function solution 



(2.15) 



(2.16) 



In the large scale limit /c — > 0, fl2.13p becomes 



SXo + 3^Xo + P (e-''^" - 1) Sxo = , 



Then, the solution is easily found to be 



Sxoin) = e-3«/2 



(2.17) 



(2.18) 



where hI'^ and H^'^ are the Hankel function of first and second kind, respectively, and are 



complex conjugate to each other, ci and C2 are constants to be determined, and 



(2.19) 



C. Large scale modes: k <^ 



Now let us consider the long wavelength modes k k^ which are already on super-horizon 
scales by the time of the waterfall transition. For these modes, we match the WKB solution 
to the Hankel function solution at some time well before the waterfall, n < and |n| ^ 1. 

In the limit k/kc — )• 0, the WKB solution (I2.15P becomes 



OYk — > — , e exp I e 



(2.20) 



where we have assumed e 



-2rn 



> 1. 
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As for the Hankel function solution, the argument is very large in the limit —n ^ 1, 
\^e~^'"'/r ^ z/e~^" ^ u. Thus using the asymptotic form of the Hankel function, 

2" 



we find that f l2.18p becomes 

-n>l 




— exp 

2>i' V vrz 



(2.21) 



-1/2 



(2.22) 

Tx \ r J \ J 

(2) 

where for notational simplicity we have omitted the term proportional to HI whose coef- 
ficient is C2. 

Comparing (I2.22p with (12.201) . we see that h\)-^ gives the correct phase factor dependence 
of (12.201) and thus we have C2 = and 



Cl 



TT if, 



gi(jy7r/2+7r/4) 



2 v/2r^ 

Thus, the long wavelength positive frequency function is given by 



fc<Cfcc 



(2.23) 



(2.24) 



2 y/2^ - V 

Here let us evaluate the mode function at the moment of waterfall n = 0. At n = 0, 
remembering that /3 ^ 1, the Hankel function takes the form 

E^)^ (^) - H'^' (-) , (2.25) 

with V ^ y/P/r. That is, the index and the argument of the Hankel function are the same. 
In this case, the Hankel function solution takes the form 
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1/3 



(2.26) 



vj 3r(2/3) 

Then, denoting by a subscript L the long wavelength modes which are on super-horizon 
scales at n = 0, we can write 

27? 



6xL{n = 0) 



exp 









il 







(2.27) 



fc^o 32/3r(2/3) y2fcfai/3 
where the numerical factor reads 2-^/?/ [3^/'^r(2/3)] ~ 1.25854, and we have defined^ 

(2.28) 



a 



Our a is equal to in Ref. 
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As we will see in the next section, we must require a ^ 1. The above result f l2.27p implies 
that all the super-horizon modes have the same amplitude at the moment of waterfall. The 
the moment of waterfall will be taken as the "initial" time to estimate the contribution of 
5x to the curvature perturbation IZc- 



D. Small scale modes: k ^ kc 

For the modes that are still on sub-horizon scales at the time of waterfall, k ^ kc, the 
WKB solution is valid until n = 0. Denoting them by a subscript S, f l2.15p readily gives 

5xs = ^-e-'^ exp (t^e-A , (2.29) 



^kc V k 

so that at the moment of waterfall 

Ho 



6xs{n = 0) = ^=-e^'/'= • (2.30) 



This is the "initial" amplitude of the sub-horizon modes. 



Before we move on, we mention that the initial amplitudes of large scale limit fl2.27p and 
that of small scale limit fl2.30p do not match at k = kc if we extrapolate from both sides, but 
are different by a suppression factor a^^^^. This indicates that in the intermediate regime 
around k = kc these two extreme values are deviating from the limiting values and smoothly 
connected jsj. In particular, this implies that the sub- horizon modes with k > kc have 
slightly different initial amplitudes from fl2.30p . However this will not affect our subsequent 
discussions because of the phase volume ~ k^ that gives rise to a sharp peak in the spectrum 
at A; ~ akc ^ kc, as we will see below. Hence we just use (12.270 for the initial amplitude of 
the large scale modes with k < kc and ( I2.30p for that of the small scale modes with k > kc- 



E. Evolution of the relevant modes after waterfall 



Having found the "initial" amplitudes of both large scale and small scale modes, now we 
can calculate the subsequent evolution of the modes until the end of inflation. 

Let us first consider the large scale modes. The solution is given by (12.240 and is valid 
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for n > as well. Then, using the asymptotic form of the large v = ^f^jr we can find 13| 

(zye-'-") = a/^cxp ( -arv"^ - - log n\ e'^^/^ _ (2.31) 
V vra \3 4 / 

This is obtained with rn <^ 1, and is thus valid for n ^ 1/r. For any sensible model of 
hybrid inflation r <^ 1, while the number of e- folds after waterfall until the end of inflation, 
Hf, is 0{1) or at most a few. Hence this asymptotic form is valid until the end of inflation. 
Then, plugging f l2.3ip into fl2.24p . we can find that after waterfall the mode function on 
super-horizon scales evolves as^ 

|5yr(n)| = — , ° exp ( -an^^'^ — -n logn | 

' ^ ^' ^2^ ^ Vs 2 4 ^ y 

= \SXLin = 0)1 ^^^^^a-^/^xp Qcn^/^ - - ^ logn) , (2.32) 

where the initial amplitude of the large scale modes dxhin = 0) is given by f l2.27p . As 
the logarithmic term indicates, f l2.32p does not hold precisely at n = but is valid for, as 
mentioned above, some time after waterfall till the end of inflation. As we will evaluate Sx 
at the end of inflation nj = C(l), we can justifiably use (12.32^ to calculate the curvature 
perturbation. 

Let us now turn to the small scale modes. An important point to calculate the evolution 
of sub-horizon modes is that the end of inflation is determined by the quanta of x which 
become tachyonic right after waterfall [7|, and affect the effective mass of in the form 
g'^i^x^). The modes which become tachyonic satisfy, by definition, [k/k^"^ < in (12.131) . 
Assuming n = 0{1), we have \f3\ ~ 2(3rn ~ a^. Hence we find that the modes with 

^ < « (2.33) 



become tachyonic. Thus a must be much greater than unity in order to have an effective 
tachyonic instability. 

To summarize, the small scale modes of our interest, which contribute to the tachyonic 
instability and control the end of inflation, are those in the interval 

kc<k<ak,. (2.34) 



^ Note that the same dependence on the number of e-folds was found from the Airy function sohitions in 
Refs. 0, Q|- But the corresponding equation solved in these references is a particular limit of the general 
equation (|2.17p and thus so does the solution, as we show here explicitly. 
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Since {k/kcY can be neglected in comparison with /3 at leading order approximation, they 
satisfy the same equation as the equation for the large scale modes, f l2.17p . Hence the 
evolution of these modes at n > is the same as that given by (12.321) . That is. 



\^Xs{n)\ = \Sxs{n 



A exp ^^an^/^ ~ \ logra^ , (2.35) 



where we have set the overall coefficient as 

2y/7r 

III. CURVATURE PERTURBATION INDUCED BY WATERFALL FIELD 

In this section, we calculate the curvature perturbation TZc by using (I2.32p and (I2.35P in 
the context of the 6N formalism. In the 6N formalism the spacetime geometry is spatially 
smoothly varying over super-horizon scales while each Hubble horizon size region is regarded 
as a homogeneous and isotropic universe. Hence we first need to smooth over the horizon 
scale Hq^, 

Sx'in) = [5x1(0) + (5x1(0))] A^exp Qan^/^ - 3n 

SxliO) + (05^X1(0)] A' exp Qan^/^ - 3n) , (3.1) 

where 5x1(0) and hence 5x^{n) is spatially varying on super-horizon scales. Note that 
we have omitted the logarithmic dependence term on n in the exponent, which is sub- 
dominant when we evaluate at n = nj = 0{1). We have also subtracted the contribution 
from the modes with k > ak since they remain stable and behave in the same way as 
the fiat Minkowski vacuum modes, in accordance with the regularization we adopted, i.e. 
{5x^in)) = at n < 0. 

With a ^ 1, from the initial amplitudes (12.271) and (12.301) we can see that the contribution 
of sub-horizon modes is much bigger than the one from super-horizon modes if the average 
is taken. At the end of infiation we have {5x^{nf)) = m? / g"^, so that using 

(^XKO))^^. (3-2) 



which follows from (I2.30p . we have 



^ = (5x'M) = [(5xi(0)) + (5x1(0))] exp (^^anf-3n; 
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:^A^ exp - 



exp I g"^/''^ - j , (3.3) 



we find 

3/2 



exp(^-an, - 3-/ J = ^.^^ • (3.4) 

Now let us rephrase the above discussion in a form more convenient for the 6N formahsm. 
With coordinate dependence exphcit, 6x^ given by (13. ip is recast as 

5x'{n, x) = Sxlin, x) + {5xl{n)) . (3.5) 

As mentioned in the first paragraph of this section, since the smoothing is done over the 
horizon scales Hq^ , there remains no spatial coordinate dependence in ((5x|)- Meanwhile, 
we do have a spatial coordinate dependence for the modes with wavelengths longer than 
Hq^ , which is what we should take care of in the context of the 5N formalism. Neglecting 
— 3n in the exponential for simplicity since a ^ 1, splitting n = n + 5n and expanding in 
terms of (5n, fl3.1l) is written as 

^xi(o) ' 

where (5x1) ^ as discussed above. 

Now we evaluate 5n at a later time, say, at the end of inflation n = nj. Here it is 
important to note that the end of inflation is controlled by the value of 5x^ at each spatial 
point, namely, 

2 

m 

6x\nf,x) = — = {6x\nf)) . (3.7) 

Analogous to the case when the value of the inflaton fleld determines the end of inflation 
hypersurface, this condition determines the end of inflation hypersurface on which the energy 
density is uniform (at leading order approximation where the contribution of the inflaton to 
the energy density is negligible). Then, using (13. 2p and (13. 4p . we flnd 

5xi(o) 



6x {n + 6n) 



((5x1(0)) A" exp (1 + 2an'/''5n + ■ ■ ■ ) , (3.6) 



1 + 



(^l + 2an/'^6r?j , (3.8) 



(^xl(o)). 

where we have truncated at linear order in 6n. Inverting this relation, we can write the 
curvature perturbation generated between the moment of phase transition and the end of 
inflation as 

jZ(^)-Sn^ ^5xi(0^ .3 9^ 

^^^^"^ 2anyM5x|(0))- ^'''^ 
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This explicitly shows that the spectrum of TZc is determined by the spectrum of 

From the result obtained previously in Sec. Ill C\ the mode function is /c-independent 
for k < kc- This implies that the power spectrum of 6xl is white: Ps^^{k) is constant 
(in the conventional terminology used in cosmology, it is blue with the spectral index of 
4: VsxLik) = ^^/(2vr^)P5x^(fc) oc A;""^ with n = 4. See (^i^ for example). Assuming that 
the spectrum has a ultraviolet cutoff at k = k^, this implies that also has the same 
white spectrum, since the convolution of two white spectra is white. Thus apart from the 
amplitude which we will calculate below, we can already conclude that V-r oc A;'^, so that the 
spectral index is strongly blue with rin^ = 4, indicating that the curvature perturbation is 
strongly suppressed on large scales. 

Before we move to the computation of the power spectrum, let us also observe that TZc 
seems to be always negative. This can be also read from (13. 8p : although 6xl may be positive 
or negative, it appears in the form of a square in (13. 8p . So irrespective of the sign of 6xl its 
contribution is always positive. Meanwhile, the left hand side of (13. 8 p is a constant. Thus, 
to compensate the positive contribution of Sxt to make the left hand side a constant, 6n is 
always negative. Also we note that the average value of TZc is not zero. 

We will consider these issues a little further later. 

Finally, before closing this section, let us discuss constraints on the model parameters. 
First we consider the condition that comes from the fact that the initial value of Sx^ must 
be smaller than the final value of it. From (13. 2 p and (13. 7p . we find 

2 247rV 127r2 

On the other hand, for this hybrid inflation model to be viable, the amplitude of the curvature 
perturbation due to the infiaton field must not exceed the observed value, V^^^ < 10~^, 

^ -'-^ ~ Uvr^A, V WMy) - (2vr)V/3 ' ^^'^'^ 



hence 



/<(27r)^10-V/3. (3.13) 



We see that both (13. lip and (I3.13P can be safely satisfied for reasonable values of the 
parameters. As a typical example, consider the case r = m^/(3ifo) = 1/10 P = 
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M'^/Hq = 100, which imphes = 20. In this case (13. lip gives 1 while we have 

(7^ < 4 X 10^^ from (I3.13p . Thus the condition (13. lip is well satisfied in this typical case. 

Let us also consider the other conditions on the model parameters. The requirement (12. 6p 
that the linear approximation to the equation of motion for x is valid implies the condition 
on A as, using (I3.13p . 

,«»!^<(^10-V^^ (3.14) 

m"^ 3 

This gives A ^ 10^^ for r = 1/10 and /3 = 100. Another condition of A comes from the 
observational constraint on the amplitude of tensor perturbations, iJ^/m|i < 10^^°. In the 
present model, since H^/nip^ = Vo/(3mpj) = M''/(12Ampi), this gives the condition 

<10-9A. (3.15) 



On the other hand, from (3 = M'^/H^ we have M^/mpj = 12A//3. Therefore we must have 

X<10-''P'. (3.16) 

Comparing with (13.140 . we see that this condition is also well satisfied for typical values of 
the model parameters. 

IV. CORRELATION FUNCTIONS 
A. Power spectrum 

In this section, we drop the subscript L from 5xl for notational simplicity. Since (TZc) 0, 
it is more relevant to consider TZc—{Tlc) rather than TZ^ itself given by (13.90 . Nevertheless, the 
difference becomes irrelevant in the Fourier space as long as we focus on a finite wavenumber. 
We will discuss this point in the next section. 

Moving to the Fourier space, we can write 

so that the power spectrum is written as 

{i^x').m,)^ (4-2) 



1 _ IGtt^ 
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where we have used (13. 2p in the last equahty. 

Before waterfall, 5x is purely quantum and it can be expressed in terms of the creation 
and annihilation operators aj^. and ak as 

= / w?'""''^" = / ik"""' ^ ■ 

where a\ and satisfy the canonical commutation relations 

[afc,at] = (277)35(3) (fc-q), (4.4) 

otherwise zero, and the mode function Xk follows the same equation as that of 5x- Since 
the Fourier component of 5x^ is written as a convolution 

= / ^^X<?5Xfc-, , (4.5) 
we have to correlate four creation and annihilation operators with different momenta, 

{{^X')k i^x')^) = 1 1^ {{SXpSx,-p) {5xi5x,-i)) . (4.6) 

To calculate the above, we should note that what we are interested in are connected 
graphs, correlating different (5x^)fc's. Thus the meaningful contractions are 

'^Xq-i) ) + ( {^Xp^Xk-p) i^Xi^Xq-i) ) , (4.7) 

while the remaining possible contractions are within the same (5x^)fc's and hence are irrel- 
evant. Then, we can easily find 

{i^x'), {5x'),) =XpXlk^plX:Xi,^iii2nf' 

X [5(3) (p + q- Z)5(3) (k-p + l) + 5(3) (p + Z)5(3) ^k-p+q-l)] . 

(4.8) 

Thus, eliminating one of the momenta using the delta functions, and using the remaining 
delta function 5^'^\k + q) to replace q with —k, we find 

{i^x'),{Sx')^) = 2 I d'p\Xp\' \x\k^p\\' S'^'\k + q) . (4.9) 

However, from (I2.27p . we have already seen that the super- horizon mode Xk is independent 
of k, and thus can be pulled out of the integral. Hence, we only have to integrate over the 
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relevant super-horizon scale momentum, for which the upper limit is k = k^.^ Therefore, 
using fl2.27p . we finally obtain 

Since this expression has, as it should, the correct delta function dependence, we can 
readily extract the power spectrum Vn- Noting (14 .2^ we find 

= - 3n/3 [r(2/3)r n, {kj ' ^''''^ 

where the numerical coefficient reads 327r^/ |3^^/^ [r(2/3)]^} ^ 1.67255. Thus, with rif = 
(9(1), the maximum amplitude is found at k = kc as V-r. ~ which is already much 

smaller than unity for a ^ 1. For larger scales, it is exponentially suppressed and thus 
becomes absolutely negligible: for example, for a scale that exited the horizon at 50 e-folds 
before waterfall, it is suppressed by a factor (e"^*^)^ ~ 10~^^. As already discussed in the 
previous section, setting Vjz oc k"''^"^, the spectrum is very blue with the index n-ji = 4. 



B. Bispectrum 

Having found the curvature perturbation and the solution of the mode function, it is 
now straightforward to calculate the three-point correlation function. We can start from the 



^ Note, however, that mathematicaUy there seems no apparent reason to set the upper hmit of the integral 
seems reasonable to extend the range of integration into sub-horizon scales up to an 
arbitrary ultraviolet cutoff at k = fcyv with /cyv 3> fcc- If proceeding with the sub-horizon mode function 
solution (|2.30p . one finds that the squared mode function |xfcP is suppressed by a factor of kc/k relative to 
super-horizon modes. But this suppression factor is not strong enough to make the integral independent 
of the ultraviolet cutoff. Since the integrand |xpP |x|fc-p||^ is proportional to p~^, the integral will be 
dominated by the contribution from the ultraviolet cutoff, leading to the result in proportional to k\jv- 
Of course there is a natural choice for the cutoff in the present case; k\jv — akc, up to which the modes 
become tachyonic, as advocated in Ref. 0, If we are to take this choice, then the resulting amplitude 
of curvature perturbations will be substantially enhanced, though the qualitative result will not change. 
Nevertheless, this strong dependence of super-horizon fluctuations on the ultraviolet cutoff deep inside 
the horizon looks physically strange because it seems to imply the violation of causality. In fact if this 
were indeed the case, then we would have a first example in which the 5N formalism fails even for the 
curvature perturbation on super-horizon scales. This may be originated from our assumption of the 
knowledge of the entire universe beyond the horizon scale in the Fourier transformation. We discuss this 
point of maintaining causality regarding the horizon scale patches in the inflating universe in Appendix iBl 
justifying (13.11) and (14. 9p which are the very foundation of our computation of the correlation functions. 
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definition 



(7^e(fel)7^e(A;2)7^c(fc3)) ={^7r fS^'\ki + k2 + k3)B^^{kl, fes, 




X ((^Xqi^Xfel-qi) (5X925Xfc2-q2) (^Xqa^Xfca 



3 -(73)) ■ 

(4.12) 



As before, we are interested in the connected graphs. This means we only take con- 
tractions between those coming from different {5x)kS. It is immediately seen that there 
are 8 possible contractions: for one of the two ^x^'s in TZdki), there are four choices of 
contractions to one of Sxk^ in 7^0(^2) and lZc{k^), and for the remaining 5xk in TZciki), 
there are two ways of contraction to either TZc{k2) or TZc{ks) which are not chosen by the 
first contraction. This gives the total number of 4 x 2 = 8 different contractions. These are 
explicitly written as 

I 1 ~i I ~l 1 I 1 ~i I ~i I 



{{^Xqi^Xk,-qi){SXqJXk,-g,)i^Xc,,5Xk3~q3)) + ((<^Xqi(^Xfel-qi ) (^Xq/^fc^-qJ (^Xqaf^Xfea-^/s)) 



+ {{^XqiSXk,-qi){SXq2^Xk2-q2){^Xq3^Xk3-q3)) + ii^Xqi^Xk^-q^) {SXqJ X k>-q,) ((hqJXk3-q3)) 



+ {{^Xqi (^Xfel -91 ) {SXq2 ^Xk2-q2 ) {^Xq3 ^Xk3-q3 ) ) + ( (<^Xqi ^Xki -91 ) {^Xq2 ^Xk2 -q2 ) {^Xq3 ^Xk3 -q3 ) ) 



+ {{^Xqi^Xki-qi){SXq2^Xk2-q2){^Xq3^Xk3-q3)) + ( (<^Xqi<^Xfel-91 ) (<^Xq2<^Xfe2-q2 ) (<^Xq3<^Xfe3-93 )) ■ 



Each of these terms exactly corresponds to the term with 6^^^ {ki + ^2 + fcs). Indeed, we find 
these 8 contractions give 



Again noting that Ixkl"^ is independent of momentum, and has a cut-off at A; = kc, we obtain 



(4.13) 




(4.14) 




(4.15) 
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Comparing this expression with the definition of the bispectrum, we find 

Bn{ki,k2,k^) = ^ , , (4.16) 

^ ^ 35[r(2/3)f fcf ^ ' 

where the numerical coefficient reads 16(27r)^/ {3^[r(2/3)]®} ~ 4128.89. To leading order, 
the bispectrum has no momentum dependence, and thus the dimensionless shape function 
{kik2k3)'^B-ji{ki, k2, ks) exhibits its maximum amplitude at the equilateral limit ki = k2 = 
k^. This is anticipated, since the curvature perturbation produced by the waterfall field 
is intrinsically highly non-Gaussian. Note, however, that this bispectrum is completely 
unobservable on large scales: in the equilateral limit, multiplying k^, we see for example 
that it is exponentially suppressed by a factor of (e"^*^)^ ~ 10"^^*^ for a scale that exited 
the horizon at 50 e-folds before the waterfall. Thus this bispectrum is totally hopeless to be 
detected on large scales. 



V. DISTRIBUTION OF CURVATURE PERTURBATION 

In this section, we consider the one-point probability distribution function of TZc- Basi- 
cally, we can guess the form of the probability distribution function. At leading order TZc 
is proportional to the square of 6xl which is very close to Gaussian. Thus, the probability 
distribution of TZc ~ 5xi is expected to be very close to the chi-squared distribution. 

The fully nonlinear distribution function can be obtained from (13. ip . By setting n = 
Hf + 6n and the left hand side of it to be m'^/g'^, and regarding 6n = TZc as a function of 
5xl, the distribution function P of IZc is given as 

P(7^,)=P,(5xL)^. (5.1) 
Here, we already know that P^(5xl) is a Gaussian distribution with zero mean, i.e. 

Px(^Xl) = ^^^exp {-^\ , (5.2) 

and the variance = ~ {^Xl)"^, with {Sxl) = 0, is given by 



2 



[5x1) = J d\ogkVsx,{k) , (5.3) 



where the power spectrum Vsxl (^) of the fluctuations 6xl can be found from f l2.27p as 

^--^^^^^'^^-'^ 3^/3.[r(2/3)f U) • ^'-'^ 
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Noting that the large scale modes has a cut-off at /c = kc, we obtain 

9 a~ 



-2/3^2 



Now, we evaluate ( 13.ip at n = + 5n = n/ + T^c to write 



(5.5) 



(5.6) 



This equation can be easily solved for 6xl as a function of Tic, 



exp 



4 



-a 



(n/ + 7^e)3/2 + 3(r^^ + 7^,) 



87r2 



(5.7) 



Thus from (15. ip we can immediately find the probability distribution of TZc as 



P(K 



2710^ 



■ exp 



XL 



:e 3^ 



X 



exp 



m 



--a(r^^ + 7^,)3/2 + 3(n/ + 7^, 
-^a(n/ + 7^e)^/' + 3(n/ + 7^e) 



. -1/2 



Stt- 



(5.8) 



This is a fairly complex probability distribution function, and is very different from the 
Gaussian one. We plot it in Fig. [TJ 

Having the distribution function of TZc at hand, let us consider the mean value (TZc). We 
can formally write it as 



(7^,) = J 7^,p(7^e)ci7^e = J 7^,(5xL)P(5xL)c^<5xL , (5.9) 

where TZc is now regarded as a function of Sxl- Although we cannot invert (15.71) to find 
TZc{Sxl) exactly, we can obtain an approximate expression by assuming \TZc\ ^ 1 as 



Tlr 



log 



m 



2tt2\ -1 



5x1 + 



2/3 



{2nf (6x 



3^1/2 \ H, 



f 



+ 



(5.10) 



where we have expanded in the limit ^ a"^ Hq / {Stt'^) . It is trivial to find that for a 
Gaussian distribution P(x), 

J x^F{x)dx = al, (5.11) 
and thus the average value of TZc is found, using (15. 5p . as 



An 



a 



-11/3 



37/3 [r(2/3)]' r^y^ 



(5.12) 
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FIG. 1: The probability distribution of \TZc\ (|5.8p . with r ~ 0.01 and a ~ 8. As \TZc\ becomes 
larger, the probability drops down extremely sharply. For comparison, we also show the chi- 
squared distribution (dotted line) with appropriate normalization. At small value of \TZc\, the two 
distribution functions behave in the same manner but they become different at larger |7^c|- 

The overall numerical factor is An/ {3^/=^ [T{2/3)f} ^ 0.527976. If we take nj = 0(1), the 
most important factor is its dependence on a: for a ^ 1, it is indeed very small. If we 
could have a of order unity, the mean value could become large. But as we have discussed 
in Sec. IIIEt this cannot be the case because of the condition for an efficient tachyonic 
instability. 

As we have mentioned before, the relevant curvature perturbation is not TZ^. itself but 
TZc — (TZc)- Therefore, although TZc is always negative, the true fluctuations from the mean 
value can become positive. Nevertheless, since the mean value (TZc) turns out to be very 
small, there is no chance to have a large positive fluctuation: the fluctuation is bounded 
from above as 



Recalling that TZc is negative of the gravitational potential, we can see that the curvature 
perturbation induced by the waterfall field repels matter around rather than attract. This 
implies that there would be no primordial black hole formation even on scales as small as the 
Hubble horizon scale, but rather bubbles of void may appear. But this is a highly qualitative 



7^c-(7^c)<|(7^c)|. 



(5.13) 
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argument and we need more explicit calculations, which we do not pursue in this note. 



VI. CONCLUSION 



In this paper, we have examined the contribution of the waterfall field x to the curvature 
perturbation. The waterfall field x can change the final curvature perturbation during the 
period between the moment of waterfall and the end of inflation, because x controls the 
physical processes during this time: the waterfall phase transition occurs as soon as the 
effective mass squared of x becomes negative, and the end of inflation is determined by the 
mean square fluctuations of the sub-horizon modes of x which became tachyonic after the 
waterfall transition. 

By solving the equation of Xi we have obtained for both the super- and sub-horizon modes 
the amplitudes at the moment of waterfall and time dependence until the end of inflation in 
terms of the number of e-folds. Using the 5N formalism, we have calculated both the power 
spectrum and bispectrum of the curvature perturbation induced by the waterfall field x- The 
power spectrum is steeply blue with n-n^ = 4, and the bispectrum exhibits the maximum 
amplitude at the equilateral limit. This indicates that the distribution of the curvature 
perturbation is intrinsically non-Gaussian, and we have presented the explicit form of the 
distribution function. On large scales, however, both the power spectrum and bispectrum 
are exponentially suppressed and totally negligible. 
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Appendix A: Evaluation of TZc from linear perturbation equation 

Here to check if our result based on the 6N formalism is consistent with the standard 
perturbation theory, we evaluate the curvature perturbation by using the linear perturbation 
equation for TZc- See also Ref. jsl for this approach. 

In linear theory, on super-horizon scales, it is known that the curvature perturbation on 
comoving slices TZc satisfies 

5P 

T^c = -H—^ , (Al) 

where 6Pc is the pressure perturbation on comoving slices. The comoving slice is defined by 
5T°j = 0. In the present case, this means 

5T^i = - (<j)di6(l) + xd^x) ~ -05i<50 = , (A2) 



where we have used the fact that x = ^Xl and Vi)y^^{k) oc k^, namely the fact that on 
super-horizon scales the contribution from the waterfall field to 5T^i is negligible compared 
to that from the infiaton field. That is, on super-horizon scales, the comoving slices are 
defined solely in terms of the infiaton as those on which the infiaton field is homogeneous. 

Therefore the contribution to the pressure perturbation 5Pc comes totally from the wa- 
terfall field 



As for p + P, we have 



5Pc = \h15x:l' + \{M^ - g^^)bxl . (A3) 



p + P = El [0'^ + Ws^)\ . (A4) 
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Let us evaluate 0'^ and {Sx's'^) to see which term dominates during waterfall. For 0'^ we 
have 



n) = (r0e)2e-2'^^ = i^o'^e"'™ • (A5) 



A'2/„\ _ /^A \2„-2rn _ zj2'_r „-2rn 

9 

,1 2 



For {Sx's ) we have at n > 1 

{Sx's\n)) ^ a\{6xUn)) ^ H',-^A' exp (^^an'/' - 3n^ . (A6) 
Hence the time dependent ratio of {5x's^) to 0'^ is written as 

^(^) = ^ ^bs:;- ( o "'^^^ - ) . (AT) 



Using f l3.7p and flA6p . at the end of inflation we have 

Therefore, for Uf > 1, (^x's^) becomes dominant toward the end of inflation. Using this 
result, we can rewrite f lA7p as 



R{n) 
Wf) 



R{n) = R{nj) ^ 6njexp 



(A9) 



3 

Let rieq = Uf — An be the time at which {Sx's"^) begins to dominate over 0'^. Since the 
growth rate of {Sx's"^) is very fast and the ratio R{nf) at the end of inflation ( 1A8I) is not 
so large, ~ 10 or so, the ((5x5^)-dominated stage appears only at the very near the end of 
inflation. An <C 1. Specifically, setting -R(neq) = 1, we find 

/any Zanjr 

Therefore, 0'^, which is almost constant in time, dominates over {Sx's'^) almost all the stage 
of the waterfall n < neq. 



With the above result in mind, we rewrite (lAip as 

dUc _ 5Pc 5xl{n) R{n) 



(All) 



dn p + P {5xl{n)) 1 + R{n) ' 

Since Sx'iin) / (Sx'si'iT-)) is time-independent, we can just replace it by that evaluated at 
n = 0. Then (lAlip can be expressed as 

dn, 6x1(0) R{n) 



dn (5x1(0)) l + R{n)- 



(A12) 
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The last factor on the right hand side is neghgible for n < Ucq and approximately equal to 
one for rigq < n < nj. Therefore, with the initial condition that TZc{0) = 0, it can be easily 
integrated to give 

With the identification that An ^ l/{2any^) as evaluated in ( lAlOp . this agrees with our 
result using the 6N formalism^. 



Appendix B: Short wavelength modes 

In this section, we justify f l3.ip and the integration of f l4.9p . and argue why we do not 
go beyond the horizon scale. We consider a scalar field (f){x) and decompose it into Fourier 
modes 4>{k). Let us call the modes with wavelengths smaller than the horizon size Lh = 
2tt/H the short wavelength modes and those larger than Lh the long wavelength modes. 
We assume the universe is inflating. 

When we decompose (f>{x), usually we assume we have the knowledge of the whole (in- 
finitely large) universe. That is, 

<Pix) = I -0^J{k) e*-^ ^ m= Jd'x<P{x)e-^''--. (Bl) 

If we divide the above into those composed of long wavelength modes and short wavelength 
modes, 

<P{x) = <P,{x) + <Ps{x) = [ -f^ me'"- + [ yf^ me''^ , (B2) 

Jk<H l^^J Jk>H l^^J 

then this will naturally induce a non-zero correlation between (j)s{x) and (psiv) even if the 
two points are separated at a distance larger than the horizon size, 

{<Ps{x)My)) ^ (B3) 

for |a3 — y I > Lh- Since each horizon size region should be causally unrelated during inflation, 
this result is acausal. This is apparently due to our assumption that we, i.e. the observers 
belonging to different regions of horizon size, have the knowledge of the whole universe. 

^ Note that if we faithfully integrate (IA12[) using (jA9l) . we can even recover the logarithmic correction 
factor ln(6n/) in (jAlOp as the leading order approximation of the integral. 
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Therefore, in stead of (lB2p . it is more reasonable to divide the field in such a way that 
4>s{x) aiid (psiy) will not be correlated if |a; — y| > Lh- To incorporate this prescription, 
we proceed as follows. We introduce two boxes of different size, a very large box where 
L = NLh with ^ 1 being a very large integer and the horizon size box Li\j. The large 
box would correspond to the present horizon size of the universe. 

We define Xs{^) for each horizon size box as 



(B4) 



where n = {ni,n2,n3) {rii are integers), Xi is the center of i-th box, and d(i){x) = 1 if a; is 
in the i-th horizon size region and zero otherwise. The long wavelength part is defined by 



\n\<N 

Thus we have the decomposition. 



(B5) 



[X 



4>l[x) + (psix) , 



(B6) 
(B7) 



This guarantees that there is no correlation of between two short wavelength modes that 
belong to two different horizon size regions: for \x ~ y \ > Lh, 



{x)4>s{y)) = o 



Now we take the square of (pix) and average over the horizon scale. We obtain 



(B8) 



X 



[X 



(B9) 



It is reasonable to assume that is independent of the region (i). Hence we may set 



s{k) . Then since = 1; obtain 



i k 

This agrees with (13. ip . 



'x) + Y.Vs{k) . (BIO) 
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Also, if we consider the sum on the short wavelength modes, 



^){P) 



p+k) =J2 



(p + k) 



(Bll) 



which appears in ( 14 .Op , it is apparent that this is non- vanishing only for | fc | > 2tt/Lh, because 
there exists no sum for < 2tt/Lh by definition. This means there will be no contribution 
from the short wavelength modes to the curvature perturbation on super-horizon scales. 
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